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Abstract 

The direct observation of extrasolar planets and their 
spectra is coming into reach with the new generation of 
ground-based near-IR interferometers, like the Very Large 
Telescope Interferometer (VLTI). The high contrast be- 
tween star and planet requires an excellent calibration of 
atmospheric distortions. Proposed techniques are the ob- 
servation of color-differential or closure phases. The dif- 
ferential phase, however, is only in a first order approxi- 
mation independent of atmospheric influences because of 
dispersion effects. This might prevent differential phase 
observations of extrasolar planets. The closure phase, on 
the other hand, is immune to atmospheric phase errors 
and is therefore a promising alternative. We have modeled 
the response of the closure phase instrument AMBER at 
the VLTI to realistic models of known extrasolar plane- 
tary systems taking into account their theoretical spectra 
as well as the geometry of the VLTI. We present a strat- 
egy to determine the geometry of the planetary system 
and the spectrum of the extrasolar planet from closure 
phase observations in a deterministic way without any a 
priori assumptions. We show that the nulls in the closure 
phase do only depend on the system geometry but not on 
the planetary or stellar spectra. Therefore, the geometry 
of the system can be determined by measuring the nulls 
in the closure phase and braking the remaining ambigu- 
ity due to the unknown system orientation by means of 
observations at different hour angles. Based on the known 
geometry, the planet spectrum can be directly synthesized 
from the closure phases. 



Top: theo. spectrum of 51Pegb. Bottom: closure phase 




Figure 1. Spectral signal of planet in the closure 
phase. Top panel: theoretical near-IR spectrum of the ex- 
trasolar planet orbiting 51 Peg displaying CO and H2 O ab- 
sorption bands (Sudarsky, Burrows & Hubeny 2003). Bot- 
tom panel: closure phases in milliradian calculated by us 
based on the theoretical spectrum of 51 Peg b shown above 
as well as on the stellar spectrum for a simulated obser- 
vation with the near-IR instrument AMBER/VLTI using 
the three telescopes UT1, UTS and UT4 (baseline length 
102m, 62m, 130m). It is evident that the spectrally re- 
solved closure phases contain a wealth of spectral informa- 
tion of the planet. 
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1. Introduction 

Since the first extrasolar planet candidate was discovered 
in 1995 orbiting the solar-like star 51 Peg (Mayor & Queloz 
1995), more than 130 extrasolar planets have been de- 
tected by radial velocity surveys (e.g. Mayor et al. 2003 
for recent discoveries). So far, the only observational infor- 
mation on the spectrum of an extrasolar planet has been 
obtained for the transiting system HD209458 (e.g. Char- 
bonnneau et al. 2002, Vidal-Madjar et al. 2004). 



In order to gain deeper insights into the physics of ex- 
trasolar planets and directly detect e.g. their atmospheres, 
observations with milliarcsec spatial resolution and at an 
unprecedented dynamical range are required (e.g. 10 -5 
in the near-IR for 51 Peg). Such observations are coming 
into reach with the new generation of ground-based near- 
IR interferometers, like the VLTI. Observations of high 
contrast systems, like extrasolar planets, require an ex- 
cellent calibration of atmospheric distortions. It has been 
proposed to observe extrasolar planets and their spectra 
through the differential phase or closure phase method 
(e.g. Quirrenbach & Mariotti 1997, Akeson & Swain 1999, 
Quirrenbach 2000, Lopez & Petrov 2000, Segransan 2001, 
Vannier et al. 2004). An overview of the mathematical 
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basis for the observations of intcrfcromctric phases of ex- 
trasolar planets has also been provided by Meisner (2001) 
and in particular on the web page referenced therein. 

Differential phases are measured between different wave- 
length bands, which, in a first order approximation, are 
subject to the same optical path length fluctuations intro- 
duced by atmospheric turbulence. The resulting differen- 
tial phase errors can therefore be eliminated in the data 
reduction to a certain degree. However, higher order ef- 
fects due to random dispersion in the atmosphere as well 
as in the delay lines (caused mainly by variations of the 
humidity of the air) can be a significant source of sys- 
tematic errors and might prevent one from reaching the 
necessary precision to observe differential phases of extra- 
solar planets. 

A more promising approach is the observation of clo- 
sure phases, which are measured for a closed chain of three 
or more telescopes. The closure phase $ c i is the sum of 
the individual single-baseline phases $^ measured on the 
baseline between telescope i and j of the array: 

$cl = $12 + $23 + $31- 

The closure phase is independent of atmospheric phase 
errors and allows therefore a self-calibration, which was 
first recognized by Jennison (1958). 

Closure phase observations from the ground will be 
possible with the AMBER (Astronomical Multi-BEam com- 
bineR) instrument at the VLTI operated by ESO (Petrov 
et al. 2003, Malbet et al. 2004) . AMBER is operating 
in the near-IR J, H, and K bands from 1.1 to 2.4/xm. It 
is currently undergoing commissioning at Cerro Paranal 
and is scheduled to be available for regular observations in 
2005. With a multi-beam combiner element and a fringe 
dispersing mode, it has the capability of combining the 
light of three telescopes (UTs or ATs), and performing 
spectrally resolved differential and closure phase observa- 
tions. There are three different spectral resolutions avail- 
able: low (R=A/AA=35), medium (R=500. . . 1000) and 
high (R=10000... 15000). 

We have modeled the response of the AMBER instru- 
ment at the VLTI for known extrasolar planetary systems 
taking into account their theoretical spectra as well as the 
geometry of the VLTI. The closure phase depends not only 
on the contrast ratio between planet and star but also on 
the (generally unknown) system geometry as well as on 
the interferometer geometry and the hour angle of obser- 
vation. In the following, we present a method to determine 
the planetary system geometry and the planet spectrum 
from closure phase observations in a deterministic way 
without any a priori assumptions. 

2. Modeling of closure phases of exoplanets 

We have written a code to simulate closure phase obser- 
vations with AMBER for known extrasolar planetary sys- 
tems, taking into account their theoretical spectra as well 
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Figure 2. Relation between the closure phase and 
the corresponding single-baseline phases. The figure 
shows simulated phase measurements for a system com- 
prised of two point sources with a 3 milliarcsec separation 
and a constant flux ratio of the components of 10~ 3 ob- 
served at zenith with the VLTI telescopes UT1, UTS and 
UT4-. The three top panels display the phases calculated 
for the single-baselines Byi, B23 and B31. (We note that 
they cannot be measured without a phase reference (e.g. 
PRIMA ) and are displayed here only for theoretical con- 
siderations.) The bottom panel shows the simulated closure 
phase for these three telescopes. All phases are in millirad. 
It is evident that the nulls in the closure phase are always 
also nulls in the single-baselines phases. In particular, ev- 
ery second null of a single-baseline phase is also a null in 
the closure phase. 

as the geometry of the VLTI. We use theoretical spec- 
tra, which have been calculated by Sudarsky, Burrows & 
Hubeny 2003 (see their paper for details). For the sepa- 
ration between star and planet the published semi-major 
axis is taken, whereas the orientation of the system is un- 
known and an arbitrary value was chosen for the purpose 
of the simulations. 

Figffl displays in the top panel a theoretical spectrum 
of the giant irradiated planet orbiting the solar-like star 
51 Peg. The bottom panel of Fig.Q] shows the modeled 
closure phase response to the planetary system 51 Peg 
in a simulated observation with AMBER and the three 
8 m unit telescopes UT1, UT3 and UT4 (baseline lengths 
102m, 62m, 130m). 

It can be seen from Fig^Jthat the spectrally resolved 
closure phases contain information on the planetary spec- 
trum. However, they also depend on the interferometer 
geometry, the hour angle of the observations, the spectra 
of the star, and on the planetary system geometry. The 
latter is generally unknown for radial velocity planets. 
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Further results of the simulations are shown in FigJ^J 
and FigOand are discussed in the following sections. 

3. Nulls in the closure phase 

We show in this section that the nulls in the closure phase 
are independent of the planet / star contrast ratio but do 
only depend on the planetary system geometry, i.e. sepa- 
ration and orientation of the system. Therefore, the mea- 
surement of the wavelength for which the closure phase is 
zero, is the first step in disentangling the planet spectrum 
from the closure phase. 

In a mathematical sense, the condition for single-base- 
line phases $^ to be zero is that the dot product between 
separation vector s (pointing from the star to the planet) 
and the corresponding projected baseline Bij equals mul- 
tiples of A/2: 

s*Bij = n • A/2 <$> tj = 0. (1) 

Thus, the nulls in the single-baseline phase depend 
only on the interferometer geometry Bij and the plane- 
tary system geometry s, but are independent of the planet 
and star spectra. (We note that the absolute interferomet- 
ric phase measured on a single-baseline is not an accessible 
quantity without phase referencing, as planned for exam- 
ple for the instrument PRIMA, which is under develop- 
ment for the VLTI.) 

We have shown analytically (Joergens & Quirrenbach 
2004) that there is also a simple condition for nulls in the 
closure phase: the closure phase is zero if the dot prod- 
uct between separation vector s and one of the involved 
baseline vectors Bij equals multiples of A: 

s * Bij = n • A |s| • \Bij\ flij = n ■ A (2) 

=► $ cl =^$ y =0 (3) 

ij 

with §ij being the angle between s and Bij. 

Comparison with Eqn.^ shows that for every second 
null in a single-baseline phase, the closure phase has also 
a null. 

This is illustrated by Fig|3 which displays the simu- 
lated single-baseline as well as closure phases for two point 
sources with a constant contrast ratio of 10 -3 for all wave- 
lengths. Considering such a 'flat spectrum' allows us to 
show the relation between the closure phase and the cor- 
responding individual phases without the complication of 
the complex spectral features of the source. It can be seen 
from the plot that the nulls in the closure phase are al- 
ways also nulls in one of the corresponding single-baseline 
phases. Furthermore, every second null of a single-baseline 
phase is also a null in the closure phase. 

This shows that the nulls in the closure phase are also 
independent of the planet / star contrast ratio and do only 
depend on the system geometry. 
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Figure 3. Closure phase time series for the planet 
system 51 Peg. Simulations for observations with the 
VLTI array UT1, UTS, UT4 for hour angles HA = -3, -2, - 
1, 0, +1 hr (from top to bottom). Vertical lines indicate the 
wavelength range covered by AMBER (the near-IR bands 
between 1.1 /im and 2.5 fim). 

4. Earth rotational synthesis 

We now proceed to describe an algorithm to determine the 
spectrum of a planet and the geometry of the star-planet 
system (i.e. the angular separation and position angle of 
the planet with respect to the star at the time of obser- 
vation) from closure phase observations, without any a 
priori knowledge of these quantities. The stellar spectrum 
can be presumed to be known, however, since this is easily 
measurable with a simple spectrograph. Determining the 
planet spectrum is therefore equivalent to measuring the 
planet / star contrast ratio. 

As described in the previous section, the nulls of the 
closure phase contain information about the system ge- 
ometry because they are related to the nulls of the single- 
baseline phases. However, it is not known a priori to which 
single baselines the individual nulls in the closure phase 
correspond. We can thus interpret Eqn.0 as follows: By 
determining a wavelength for which the closure phase is 
zero, we have a measure for the projection \s\ ■ #y of the 
separation vector s onto one of the individual baselines 
Bij but we do not know onto which one. Furthermore, we 
do not necessarily know the order n of the null. The geo- 
metrical locus of all possibilities for s is therefore a set of 
straight lines, as shown in the left panel of Fig.0] Each line 
in this figure is perpendicular to one of the baselines and 



Figure 4- Possible geometries of a star planet system as derived from nulls of the closure phase. Left: 
Determination of the nulls in the closure phase signal of a planetary system at a certain hour angle is a measure for 
the projection of the separation vector s between star and planet onto one of the three projected baselines. Due to the 
unknown orientation we do not know onto which one. The dashed lines indicate the possible locations of s for the 
nulls corresponding to n = 1 and 2. Right: Observations at a later hour angle gives another set of lines. The possible 
locations for s are now reduced to several grid points as marked by the intersections with the former lines. With a third 
observation, it will be possible to select the correct intersection point. 



corresponds to an assumed set of i,j,n. For clarity only 
the lines corresponding to n — 1 and n — 2 are shown. 

The set of projected baseline vectors formed by the 
three telescopes changes with time due to Earth's rota- 
tion. Therefore, a second observation gives an indepen- 
dent constraint on the geometry as shown in the right 
panel of Fig.Q] The possible system geometries are now 
the intersections of the set of lines corresponding to the 
first observation with the set corresponding to the second. 
This discrete set of points has been marked with filled cir- 
cles in the figure. A third observation will produce yet an- 
other independent set of lines. In general, this set will pass 
through only one of the previously marked points: this is 
the true separation vector. It is thus possible to derive the 
system geometry unambiguously with three observations. 

Once the separation vector is known, it is straightfor- 
ward to determine the planet / star contrast ratio for each 
wavelength through a numerical inversion of the expres- 
sion for the closure phase. The only difficulty occurs at or 
very close to the nulls, where the signal-to-noise is zero 
or very low. However, since three observations are needed 
in any case for the determination of s, one can perform 
the three inversions and combine them with (wavelength- 
dependent) weights appropriate for their respective signal- 
to-noise ratios. Consider, for example, the wavelength range 
around A = 1.67 /xm (1/A = 0.6 /im -1 ) in Fig. EJ The 
observation at — lhr (third panel) gives very low signal- 
to-noise close to the null near 1/A = 0.6 /an -1 , but ob- 
servations at other times can be used to infer the planet 
spectrum in this wavelength region. 



5. Outlook 

We have shown that the separation vector and spectrum 
of extrasolar planets can be determined from closure phase 
measurements in a deterministic way, with a non-iterative 
algorithm, and without any a priori assumptions. For a 
practical application of this technique, a few additional 
complications will have to be considered: 

— The assumption that the star and the planet are point 
sources will have to be relaxed. Taking into account 
that the star is slightly resolved by the interferome- 
ter complicates the analysis and implies that the rela- 
tion between nulls in the closure phase and the single- 
baseline phases is only approximately fulfilled. 

— For observations from the ground, the useful wave- 
length range in the near-infrared is limited to the at- 
mospheric windows and thus non-contiguous. 

— The finite signal-to-noise ratio of realistic observations 
will lead to an uncertainty in determining the exact 
wavelengths of closure phase nulls. The lines and in- 
tersection points of Fig. 0] will thus be broadened. 

— The non-zero time needed to accumulate sufficient signal- 
to-noise on the closure phase means that the interfer- 
ometer geometry will change slightly during the obser- 
vations. This is again equivalent to a slight broadening 
of the allowed regions in Fig. 01 

— If the time required to accumulate all observations is 
not short compared to the orbital period of the planet, 
the motion of the planet will also contribute an uncer- 
tainty in the derived geometry. 

These practical complications will certainly be at least 
partially compensated by a larger number of observations; 
one would probably take ten or more rather than the min- 
imum three. One could also combine the algorithm pre- 
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sented here with a global x 2 minimization, in which the 
orbital parameters and spectrum of the planet are fitted 
to the observations. The purpose of our algorithm would 
then consist of providing a robust starting point for the x 2 
minimization algorithm, which is usually crucial to ensure 
convergence to the correct minimum. 

The AMBER instrument has arrived on Cerro Paranal 
this year and the first commissioning runs have taken 
place. It is now necessary to determine if the required clo- 
sure phase precision (better than 0.1 mrad, see Fig.^l can 
be reached. If so, the VLTI will provide unprecedented op- 
portunities for observations of extrasolar planets and their 
spectra. 
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